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The idea that in dynamical wave function collapse models the wave func- 
tion is superfluous is investigated. Evidence is presented for the conjecture 
that, in a model of a field theory on a 1+1 lightcone lattice, knowing the 
field configuration on the lattice back to some time in the past, allows the 
wave function or quantum state at the present moment to be calculated, to 
arbitrary accuracy so long as enough of the past field configuration is known. 



1 INTRODUCTION 

The question of the status of the state vector, Psi, in standard text-book 
quantum mechanics has been a controversial issue since Bohr's and Einstein's 
time. Is the state vector a complete description of the physical state of a 
system, or is it incomplete and needing of completion with extra informa- 
tion, or does it represent a state of (someone or something's) knowledge of 
the system? The class of models known variously as "dynamical collapse 
models" and "spontaneous localisation models" are observer independent al- 
ternatives to standard quantum theory. We can ask: what happens to this 
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question about the status of the state vector? Is it immediately resolved by 
the formalism of the collapse models or is there still a question to answer? 

The general structure of all dynamical collapse models is similar: there is a 
state vector, Psi, which undergoes a stochastic evolution in Hilbert space and 
there is a "classical" (c-number) entity - let's call it "q-bar" following Diosi 
- with a stochastic evolution in spacetime. The stochastic dynamics for the 
two entities - Psi and q-bar - are coupled together. The stochastic dynamics 
in Hilbert space depends on which q-bar is realised in such a way as to tend 
to drive Psi into an eigenstate of an operator (q-hat) that corresponds to 
q-bar: this is the eponymous "collapse" in these models. And the probability 
distribution for the realised values of q-bar depends on Psi. 

The choice of q-bar varies from model to model. In the original GRW 
model j2] and a proposed relativistic version ^3 , q-bar is a sequence of discrete 
"collapse centres" or spacetime events, in Diosi's model for single particle 
quantum mechanics |1_ q-bar is a particle position (see, however, footnote 1), 
in Continuous Spontaneous Localisation (CSL) models |^ [S] q-bar is a scalar 
field. In all cases the c-number entity q-bar is defined on spacetime and is 
therefore covariant in essence. 

The Bell ontology [S] for the GRW model states that the collapse cen- 
tres are the beables or real variables. The analogous ontology for collapse 
models in general is that the history of q-bar - whatever it happens to be 
in the model - is real. Work by Diosi shows that any prediction about re- 
sults of macroscopic experiments and observations that can be made using 
the expectation value of operator q-hat in state Psi, can also be made. For 
All Practical Purposes (FAPP), using only knowledge about q-bar, suitably 
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regularised and coarse grained. Indeed, in non-relativistic theories q-bar is 
equal to this expectation value plus white noise with zero mean^ EH- Put 
another way, suppose one has run one's computer simulations of the collapse 
model up to a time well to the future of anything one is interested in and 
in the computer memory is a history for Psi traced out in Hilbert space and 
a (regularised) history for q-bar traced out in spacetime. If the computer 
has a memory failure and loses all information about Psi, the information 
about q-bar would be enough, when suitably coarse grained, to make all the 
macroscopic predictions that could be made from Psi. 

An example is the lattice field theory _8; that is the subject of this paper. 
In Pj it was argued that a coarse graining of q-bar - in this case q-bar is a 
{0, l}-valued field on the lattice - displays the same structure, FAPP, as the 
coarse grained expectation value of the field operator in the quantum state 
Psi. 

In taking this point of view, that q-bar is real, we are forced to address 

the question of the status of Psi. Diosi takes the view that both Psi and q-bar 

are real |10j . In this paper we will investigate the possibility, raised explicitly 

by Kent jllj for the GRW model, that Psi doesn't exist at all - that it is 

at most a convenience and conveys no information that is not carried by the 

history of q-bar itself. 

^This raises the objection that the q-bar history is not really properly defined at all as it 
contains a white noise term. One could fall back on the argument that spacetime is widely 
expected to be fundamentally discrete and this discreteness would provide a physical cutoff 
for the frequency of the white noise. Or turn the argument around and say that if the 
Bell ontology for collapse models is desirable, this suggests the necessity of fundamental 
discreteness. 
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One can argue that there is already a way partially to demote the quantum 
state in collapse models from its status as a really existing thing to that of a 
"dynamical law". This view can be taken in formalisms in which spacetime 
histories of the system are primary (including for example consistent histories 
|12l 1131 rm I15j . Sorkin's quantum measure theory jl6l I17j and Bohmian Me- 
chanics, see e.g. ^H])- If we consider a collapse model to be a stochastic law 
for the q-bar histories then the quantum state Psi can be formally relegated 
to the initial surface from which it need never evolve. The initial state gives 
us the dynamical law for the future q-bar histories in the form of the proba- 
bility distribution on them and is not itself real. However, we can, if we know 
the q-bar history up to some spacelike surface, define an "effective" quantum 
state on that spacelike surface which tells us how to calculate the probability 
distribution on q-bar events to the future of the surface conditional on the 
known past history. The "evolution" of this effective quantum state from 
surface to surface (which is precisely the stochastic process in Hilbert space 
mentioned above) is akin to a "Bayesian" updating - on the actualisation of 
stochastic events - of the rule which gives the future probability distribution 
and is not the evolution of something physical. On this view, the quantum 
state is something we invent in order to render the dynamics Markovian. 

It would be desirable to go further than this. The initial state on the initial 
surface hangs around like the smile of the Cheshire Cat - rather insubstantial 
but still persistently there. Moreover, in the quest to make a relativistic 
collapse model, the need to begin with a state defined on an initial surface 
breaks Lorentz invariance. In this paper we elaborate on the conjecture made 
in that in collapse models even the initial state can be eliminated as a 
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necessary part of the theory (and the only information that remains from the 
state is a classical distribution over superselection sectors). We suggest that 
it can be replaced by an "initial period of q-bar history" . Knowing this initial 
period of history would allow the calculation, FAPP, of an effective quantum 
state which could be used to make predictions from then on. 

In section 2 we briefly describe a collapse model for a field theory on a 
1+1 null lattice that we will use as a testing ground for our conjecture. In 
this model, q-bar is a field configuration of O's and I's on the lattice. In 
section 3 we state the conjecture and in section 4 we describe the simulations. 
The results reported in section 5 suggest that if the field configuration is 
known to a certain depth in time Tconverge, the state vector can be deduced 
FAPP from that configuration. Thus the evolution of the field alone would 
be approximately Markovian on time scales larger than Tconverge- Section 6 
contains a summary and discussion. 



2 CAUSAL COLLAPSE MODEL ON A LIGHT- 
CONE LATTICE 

We briefly review the spontaneous collapse model that we will use to 
investigate the conjecture. We follow the presentation of 9 and refer to that 
paper for further details. The model is a modification of a unitary QFT 
on a 1+1 null lattice, making it into a collapse model by introducing local 
"hits" driving the state into field eigenstates. The spacetime lattice is N 
vertices wide and periodic in space, extends to the infinite future, and the 
links between the lattice points are left or right going null rays. A spacelike 
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surface a is specified by a sequence of N leftgoing links and N rightgoing 
links cut by the surface; examples of spatial surfaces are shown in figure 1. 
We assume an initial spacelike surface uq. 

An assignment of labels to the vertices to the future of ctq, ^1,^2, ... , is 
called "natural" \i i < j whenever vi is to the causal past of Vj. A natural 
labelling is equivalent to a linear extension of the (partial) causal order of 
the vertices. A natural labelling, vi,V2-, - ■ ■ is also equivalent to a sequence of 
spatial surfaces, ai, (T2, . . . where the surface Uk is defined such that between it 
and ctQ) lie exactly the vertices vi,. . .Vk- One can think of the natural labelling 
as giving an "evolution" rule for the spacelike surfaces: as each vertex event 
Vk occurs, the surface creeps forward by one "elementary motion" across 
that vertex. For any natural labelling and any fc, the finite set of vertices 
{vi,V2, ■ ■ ■ Vk\ is a stem, a finite set that contains its own causal past. 

The local field variables a live on the links. These field variables take only 
two values {0, 1}, so that on each link there is a qubit Hilbert space spanned 
by these two states. We denote by {aR^^ai^} (atj^ for short) the values of 
the field variables on the two outgoing links (to the right R and to the left L) 
from vertex Vk- (Note, in paper 1 we used the hatted symbol a to denote the 
actual value of the field variable but here we will use the unhatted a.) One 
can, colloquially, consider the field values and 1 to represent the absence or 
presence (resp.) of "bare particles" on the lattice. 

A quantum state \ipn) on surface is an element of the 2^^ dimensional 
Hilbert space which is a tensor product of the 2N 2-dimensional Hilbert 
spaces on each link cut by cr„. The basis vectors (the "preferred basis") of this 
Hilbert space are labelled by the possible field configurations on (T„, namely 
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the 2A^-element bit strings {0, 1}^ . We will often refer to the number of I's in 
the bit string labelling an eigenstate as the number of particles in that state. 
The Hilbert space is the direct sum of 2N + 1 sectors each of fixed particle 
number. We identify the Hilbert spaces on different surfaces in the obvious 
way using the field basis. At each vertex Vk, there is a local evolution law 
which is given by a 4-dimensional unitary "R- matrix" U{vk) (4-dimensional 
because it evolves from the two ingoing links the two outgoing links). For 
this paper we choose these R-matrices to be uniform across the lattice. (One 
can simulate external interventions by fiddling with the R-matrices.) 

In the standard text-book unitary theory, one postulates the existence of 
an external measuring agent and then this formalism can be used to predict 
the results of sequences of measurements of the field. One way to do this is 
to identify projectors: P(Qt,j.) projects onto the subspace of the Hilbert space 
spanned by the basis vectors in which the field values at vertex Vk are a^^. 
(recall there are two links outgoing from Vk and so ay^. is really two values). 

Then the joint probability distribution for the agent to measure a partic- 
ular field configuration {a^,^, a^j) • • • o^d^} on the lattice between the hyper- 
surfaces (To and an is: 

pstandard QM(^^^^^^^^ _ _ _ ^ \\P{a,jU{vn) . . . P{a,,)U{v,)\i;o)\f . 

(1) 

This probability rule is independent of the linear ordering of the vertices 
and depends only on their causal order. The rule evades the potential danger 
of violating relativistic causality described in in two ways: the causal 
structure of the vertices of the lattice is a partial order from the start (no 
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transitive completion is required) and the field variables being measured are 
completely local quantities. 

Inspired by the GRW model with Bell's ontology, this unitary quantum 
field theory requiring external agents can be turned into an observer inde- 
pendent theory for a closed system by replacing the projection operators for 
measurements in by positive operators (for "unsharp measurements") and 
adopting the resulting formula as the probability that the corresponding field 
configuration occurs. More precisely, we define on each link (i.e. on each 
2-dimensional Hilbert space associated with a link) the two operators Jq and 
Ji where 

with < X < 1. Note that Jq + Jf = 1 and this is a positive operator valued 
measure. Then we define the jump operator J{ayi^) on the 4-dimensional 
Hilbert space on the outgoing links from as the tensor product of the two 
relevant 2-dimensional jump operators, e.g. when 0^,^. = (0,0), J{a^^) = 
Jq (8> Jo- We promote Jia^^) to an operator on the Hilbert space of any 
spatial surface containing those two links by taking the tensor product with 
the identity operators on all the other components of the full Hilbert space. 
The probability of the field configuration {a^j, . . . a^^} is given by 

P(a,,,...a,J = \\J{a,jU{vn)...J{a,,)U{vi)\^o)f • (3) 
From this we can see the importance of the fact the jump operators form a 
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positive operator valued measure, which ensures consistency: 

P(a^,i,...a^,„_J = ^P{ay,,...ayJ (4) 

= ^{M---Ji(^v„)JiavJ-..\^o) (5) 

= \\J{a,„.,)U{v^-i)...Jia,,)U{v,Mo)f • (6) 

Again, (jSJ depends only on the (partial) causal order of the vertices be- 
cause any other choice of natural labelling of the same vertices gives the same 
result. These probabilities of the field configurations on all stems are enough, 
via the standard methods of measure theory, to define a unique probability 
measure on the sample space of all field configurations on the semi-infinite 
lattice. 

We stress that whereas equation Q is the probability for measuring a 
particular field configuration in standard unitary quantum theory, equation 
(jSJ is interpreted, in the Bell ontology, as the probability for the field to 
be in that configuration. The full content of the theory is the probability 
distribution © on possible field configurations, dependent on an initial state 
IV'o). 

The state on the hypersurface cr^ that is reached after the elementary 
motions over vertices vi, . . .Vn and the field values {0^,^, . . . have been 
realised is the normalised state 

II y J{(^vr.)U{vn) ■ ■ ■ J{av^)U{vi)\'^Q) . . 

'^"^ ||J(a,jC/K)---^Ki)f^(«i)lV'o)ir ^ ^ 

Thus, the probability for state (jT)) on hypersurface (T„ is ©. In order to make 
predictions about the field on the lattice to the future of - conditional on 
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the past values {ck^^ , . . . } - it is sufficient to know \ipn)- Indeed, tfie 
conditional probability of {a^„_,_i, . . . civ„+m} given by 

) = \\J{avn+m)U{Vn+m ■ ■ ■ J{avr,+i)U{Vn+l)\ll^n)f ■ (8) 

The state vector provides these conditional probabilities, and is therefore a 
convenient way of keeping the probability distribution up to date, given past 
events. 

3 THE STATUS OF THE WAVE FUNCTION: "AN 
EXECUTIVE SUMMARY"? 

We are interested in investigating the possibility of doing away with the quan- 
tum state entirely as a fundamental concept in collapse models and will be 
using the lattice model described above as a test case. As mentioned in the 
introduction, we can relegate the quantum state to a state, {tpo), on an initial 
surface ctq from where it acts as a "dynamical law" , specifying the probability 
distribution on field configurations to the future of ctq. Can we weaken even 
this status? 

We make the conjecture that, if the field configuration is known between ctq 
and an, then even if the state on ctq is not known, the state on (7„ is calculable 
up to a correction that goes to zero as n — oo. This would mean that although 
the evolved state on (j„ a priori depends on both the initial state on ctq and 
on the field values that have actually occurred in between, its dependence on 
[■^o) dies away as time goes on until all we need to know to make predictions, 
FAPP, is the field configuration back to a certain depth in time. We would 
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then have an interpretation not only assigning reahty to a field configuration 
in spacetime but further demoting the wave function by denying it a role as a 
necessary entity to the theory: the state \ipn) can be deduced FAPP from the 
field configuration to the past of (t„ to some depth in time (or exactly if the 
whole infinite past history is known) and becomes an "executive summary" 
of the past reality containing no independent information. 

More precisely, let |^q) and be two states on (Jq. Then they give, 
according to two probability distributions. Pi and P2 on field config- 
urations to the future of fJo- Choose any linear ordering of the vertices to 
the future of do, fi, ^2, • • • • Adopt the notation a(n) for a field configuration 
between ao and (T„, and |V'^^,a(n)) for the state on an that arises from jV'o) 
on (To after a(n) has happened (a = 1, 2). 

Conjecture: There exists a complex phase A such that 

|||V'n,a('^)> - AlV'n; «("-)) II ^ as n ^ 00 (9) 

for all a{n) except those which almost surely do not occur according to both 
Pi and P2. 

Conjecture (density matrix form): 

II Pi(a(n))|Vi,a(n))(V'n,a(n)| 

a{n) 

- XI Pi(«("))l^n>"("')>(V'^)"(")lll ^ as n ^ 00 (10) 

a{n) 

where || • || is the operator norm, and similarly with 1 and 2 interchanged. 
Note that we already know that the conjectures cannot be true strictly 
as stated because of the possible existence of "superselection sectors" in the 
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Hilbert space. For example, the jump operators J preserve particle number 
and if the R-matriccs do so also (this is the case we will study in detail in 
the next section) then a state in the fc-particle sector can never approach a 
state in the Z-particle sector if k ^ I. If the R- matrices preserved only particle 
number mod-2 (by allowing pair creation and annihilation of particles) then 
there would be two superselection sectors (even and odd particle number) . It 
should be noted that even if there is a conserved quantity - particle number, 
say - this quantity is conserved in the state vector but not in the realised 
field configuration. We expect however that the "conservation law" will be 
reflected in the probability measure in the sense that a suitable property of 
the coarse grained field configuration will be predicted with probability close 
to one. 

When there are superselection sectors, an initial quantum state corre- 
sponds to a classical probability distribution over the sectors and a quantum 
state in each sector, in the familiar way. Without loss of generality therefore 
we will assume in what follows that we are restricted to a single superselection 
sector and the conjectures apply to each superselection sector individually be- 
coming, effectively: two states in the same superselection sector tend to each 
other up to a phase for all histories except for a set of histories which has 
measure zero in the probability measure of both states. 
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4 THE SIMULATIONS 



We sought evidence for the conjecture in the following way. We chose a 
unitary R-matrix, uniform across the lattice, of the following form: 
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This gives a particle number preserving dynamics for the state, since the hit 
operators J also preserve particle number. 

We chose ctq to be a constant time surface and we chose two initial states, 
I^q) and \iPq) (in the same superselection sector, which here meant the same 
particle number sector). We generated, at random according to the probabil- 
ity distribution Pi or P2 field configurations to the future of do. For each of 
these field configurations, a(M) (where M was large enough for the calcula- 
tion in hand) we calculated the two states a{n)), a = 1, 2 on the surface 
an which is the n*'^ surface in a sequence of surfaces chosen according to the 
stochastic rule "choose the next elementary motion at random with uniform 
probability from those possible." This is not a covariant rule - it is equivalent 
to a probability distribution on linear extensions of the partial order on the 
whole future lattice but it does not give each equal weight - and moreover 
a covariant, Markovian rule does exist 20^ but we made the choice for ease 
of calculation. We will comment on what significance this has for our results 
below. 
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We would like to show that the two states \^p^,a{n)), a = 1,2 become 
close, up to a phase, as n gets large and more precisely we would like to know 
how the difference behaves with n. One can argue that it is not the states 
themselves that should be compared but the probability distributions for the 
field variables that they produce. Indeed, in an interpretation in which only 
the field is real, it is only this probability distribution and not the state itself 
which has physical import. 

In principle the entities that should be compared are the two probability 
distributions, for states 1 and 2, over field configurations to the future of 
any surface a, given the values of a(M) lying to the past of a. This is 
calculationally impractical and we used two simplifying strategies. First, we 
sampled the space of all surfaces by choosing a sequence ai,a2, ■ ■ ■ according 
to the rule described above. This rule does not sample uniformly in the space 
of surfaces, as mentioned above, and an improvement of our scheme would be 
to determine and then implement the covariant rule which does. Second, we 
compared the probabilities, not for the whole future field configuration but 
only for the value 1 on each of the two outgoing links from Vn for each n. 

An important point is that, as discussed in jOj, the interesting physical 
regime for these models is when the parameter X is close to one, alternatively 
when e = 1 — X is close to zero. This means that the hits are very gentle and 
superpositions of microscopically different states will last for a long time. In 
this case, however, the conditional probability of a 1 on each link becomes 
very close to 1/2, indeed it is equal to 1/2 + 0(e). (Here we clearly see the 
white noise term in the field configuration that is to be expected from Diosi's 
work.) So for small epsilon the probabilities will be close, whether or not 
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the states are coming close to each other. Indeed, let the link in question be 
denoted / and suppose, at some stage in the dynamics, / is one of the outgoing 
links from the vertex that has just been evolved over. Let the state on the 
current spacelike surface through / be denoted schematically by 

1^) =a|0) + 6|l) (12) 

where |0) (|1)) is short hand for the normalised superposition of all the terms 
in the state in which the value of the field on Hs (1). The probability that 
the field will be 1 on / (conditional on the past evolution to that stage) is 



1 + X2 1 + X2 ' ' 1 + X2 ■ 

So, for the difference between the probability for a 1 on link / in state 1 and 
in state 2 we will obtain: 

{\hi?-\h2\')\^,. (14) 
When e is small this becomes 

m^-\h2?){e + 0{e^)). (15) 

Prom this we see that the appropriate quantity to calculate for each link 
is l^ip — |62p: that gives a measure of the difference of the probability distri- 
butions that affects the coarse grained, renormalised field configuration (see 
[S]) and indeed it is a measure of the difference between the states themselves. 

Thus, we calculated for each vertex Vn (recall the sequence vi,V2,... is 
equivalent to a linear extension and is the one chosen at random by our 
evolution rule described above) and for each outgoing link, /, from Vn, the 
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quantity \bi\^ — which we denote by B{1). Of course this is a very 
approximate measure of the difference between the two states and to overcome 
this, we took the sum of this quantity over every hnk in a block, a certain 
number, m, of lattice time steps long and the width of the whole spatial 
lattice: 



where the sum is over all links with lattice time coordinate from t through 
t + m—1 (the lattice time step is 1). Our convergence criterion was Bjnit) < 5 
and we define the convergence time Tc to be the smallest time such that 



With the help of numerical simulations on 8, 9 and 10 vertex lattices we 
studied the dependence of Tc on e, on particle number, on 9 and on different 
types of initial state within fixed particle number sectors. In total about 
600 simulations were run. We also studied the convergence of states for field 
configurations not generated according to the probability distributions from 
either state, for example the field configuration (a) of all I's, (b) of all O's 
and (c) randomly generated with uniform probability distribution of 1/2 for 
a 1 on each link. We failed to find convergence only in the cases (a) and (b) 
mentioned above when the field configuration was all I's or all O's which is 
consistent with the conjecture because they almost surely do not occur in Pi 
and in P2 . Convergence occurred but was slower for the field configurations of 
type (c) than for those generated by (and therefore likely in) the probability 
distibutions of states 1 or 2. 

In our simulations we were limited as to lattice size by the exponential 
growth of the problem in vertex number, and it is at present unclear whether 




(16) 



Bruit) <6,yt> Te. 
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the limited size of the lattice has important implications for our results, in 
particular the question whether the periodic boundary conditions of the lattice 
stimulate convergence remains open. 

5 Results 

To begin by giving a flavour of the kind of simulations run, the convergence 
of two initial states is illustrated in figure [21 Each cell corresponds to a single 
link of the lattice (so there are twice as many cells across the lattice width as 
vertices) and the darkness of the cell is (positively) proportional to |6p (see 
equation p2() '). In plot (a) we show the evolution of state 1, which begins 
as an eigenstate with 4 particles on the left hand side of the lattice in the 
leftmost panel, time proceeds up the page and then the lattice continues at 
the bottom of the next panel and so on. Plot (b) shows the evolution of state 
2, which begins as a state with 4 particles on the right hand side of the lattice. 
The parameters for the evolution are X = 0.65 and 6 = 0.267r. Comparing 
(a) and (b) it can be seen that the plots are indistinguishable from halfway 
up the first panel in each. The plots after this time are somewhat superfluous 
but we show them to emphasize that we checked that the convergence persists 
long after our convergence criterion is reached. 

FigureElshows the quantity i?io(t) defined in equation (|16|) plotted against 
lattice time t, for an 8 vertex run with X = 0.95, ^ = O.lvr and initial states 
which are two different one-particle eigenstates. It shows a pleasingly sharp 
falloff to zero. 

Figure|l]is a plot of log(Tc) against log(e) for many 8 vertex runs of varying 
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e, where we chose -Bio ^-s our measure of difference and 6 = 10~^ to define 
the convergence time. The other parameters were 9 = O.lvr and {ipi) and 
\ip2) are two fixed one-particle eigenstates. The field configuration is chosen 
according to the probability distribution from \tpi). The plot is consistent 
with a dependence of 

oc ^ (17) 

As a consequence of this, in a continuum limit in which the lattice spacing 
a ^ and e = 0{^/a), the "physical" convergence time, aTc would tend to 
some finite non-zero value. 

A plot of the convergence time (defined by Biq and 5 = 10"'^) against 6 
is given in figure 121 for X = 0.95 and fixed initial one-particle eigenstates. 

The plot is difficult to interpret. It seems particularly odd when one 
realises that for 9 = it/2 and ^ = (the two limits of the range of 9 shown) 
the R-matrix is such that it does not introduce any superpositions into the 
states. Indeed the evolution is completely deterministic: for 9 = it/2 an initial 
one-particle field eigenstate remains essentially constant - just acquiring a 
phase of i at each lattice time step ~ and for 9 = it propagates at the speed 
of light along the null direction it starts off in. The evolution in these cases, 
therefore, does not "mix" the Hilbert space and states 1 and 2 can never 
converge. 

The plot, however, suggests that for values of 9 close to these limiting ones, 
convergence is faster than for 9 in the middle of the range, and indeed the 
convergence time is tending to zero. We speculate that this has something 
to do with the competing effects of the mixing by the R-matrices and the 
converging effect of the hits. If we imagine starting with two states which 
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have support over the whole of the one-particle sector of Hilbert space, the 
harder the hits, the faster the states will converge. In the extreme case, if 
e = 1 then the hit operators are projectors, state 1 will collapse into one of the 
eigenstates after one time step, the field configuration will be the one given 
by that eigenstate and state 2 will be forced into that state also. When e and 
6 vary, there is a competition between the driving towards eigenstates by the 
hits and the mixing (introduction of superpositions) by the R-matrices. In 
the runs plotted here, we kept e fixed so the strength of the hits does not vary 
but as 9 tends towards the two limiting values it could be that the R-matrix 
evolution loses the competition. If the hits drive state 1 very quickly into an 
eigenstate, then as long as there's been enough mixing so that there is even 
a tiny amplitude for that eigenstate in state 2, there will be convergence. 

Figure IHl shows results from runs on an 8 vertex lattice with X = 0.9, 
6 = 7r/4. In any given run the two initial states are eigenstates with the same 
particle number, which varies across the runs. The plot is of Tc (defined by 

and 6 = 10~^) against particle number. 

The plot is consistent with expected behaviour. The fixed particle number, 
m, sectors have dimension (2A^)!/m!(2A^ — m)! which increases as m increases 
to N and then decreases symmetrically as m increase further to 2N. When 
the Hilbert space is larger, we expect that convergence will take longer as it 
takes longer for each state to mix and acquire amplitudes for all the different 
eigenstates. We expect the plot to be symmetric because, further, there is a 
duality in the models between field value and field value 1. 
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We checked our results by taking several runs and calculating the quantity 

Cn = l-\\{il^la{n)\i^la[n))f (18) 

and comparing it to the quantity i?io(i) for the same run (recall that the way 
BiQ{t) is defined, there is one value for every tenth lattice time coordinate, 
so there are 80 times as many C„ data as i?io data). Figure is a C„ plot of 
the run shown in|21 This is a more direct check of the conjecture and in the 
future we would want to redo our analysis using this method. 

However, we present more evidence in figures 1^1 and ^2 that indicates 
that the results will be the same. Indeed, even in the details of how the 
convergence occurs in each run, the behaviours of the measures BiQ{t) and 
Cn match each other very well. On noting that the number of elementary 
motions is 8 times the lattice time, it can be seen that the main features of 
the two types of plot are well matched in time. Figures El IH show a plot of 
i?io(t) and C„ data for a run with the same initial states and parameters as 
for the simulation whose data is shown in figures |21 and [3 while figures 1101 
111! show a plot of -Bio(i) and C„ data for a run with the same states and 
parameters as for |H1 and |^ but a different = 0.257r. 

6 DISCUSSION 

We can state the import of the conjecture we have made thus: given some par- 
ticular field history from f = — ootot = then there is a physical probability 
distribution on the field histories for t > which we can express conveniently 
in the form using a quantum state at t = which is precisely specified 
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by the past field history. If one could discover an algorithm for transforming 
the data in the field history directly into the probability distribution then one 
would have built a model in terms only of the field variables which makes no 
reference to a quantum state. 

We have presented evidence for this conjecture. The falloff seen in figures 
01 and [TBTTI suggest the stronger conjecture that there is a time scale Tc such 
that even if we know only the field history from t = —Tc to t = then we 
can construct a quantum state that gives the correct predictions FAPP. We 
would like to study this further by investigating the dependence of Tc on the 
degree of convergence 6. 

For the purposes of this paper we chose 6 = 10~^ to define Tc and we 
presented evidence that Tc is of order as e ^ 0. In a continuum limit 
where the lattice spacing a — > and e = 0{y/a) then the physical convergence 
timescale, aTc would remain finite and the dynamics would be approximately 
Markovian for time scales larger than this. 

It would be valuable to check all our results by redoing the simulations 
and calculating, instead of Bmit), Cn on each sampled surface and examining 
how it tends to as we did for some runs described in the last section. 
Improvements on our methods would include calculating and implementing 
the covariant evolution rule for surfaces which would make our sampling of 
surfaces uniform. We would like to analyse quantitatively the dependence 
of Tc on the dimension of the particle number sector implied by the results 
shown in figure IHl 

Results with different types of R-matrices as well as general initial states 
are still to be investigated. In particular, further evidence for the conjecture 
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can be obtained by choosing pair-particle conserving matrices, as well as 
general matrices with no conservation laws. 

We stress that the analysis and simulations presented in this paper are at 
a rather mathematical level. The question of physics has not been addressed. 
This would involve the settling of the issue of the competition between the 
R-matrices and the hits in the collapse of superpositions of eigenstates 0. 
This bears on the conclusions of the current paper. The physically interest- 
ing range of parameters is when e is very small and 9 is also small so that 
"microscopic" superpositions persist for a long while but eventually collapse. 
In this regime, the hits are very gentle and the "mixing" of the Hilbert space 
by the R-matrices is slow. Investigating this regime is essential if we are to 
draw physically relevant conclusions about collapse models of this sort. 

Finally we extend our conjecture to all collapse models. It would be 
interesting to study it in other cases such as the GRW model and Diosi's 
single particle model. 
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Figure 1: The light cone lattice, at is a constant time surface and cr is a 
generic spacelike surface. 
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Figure 2: (a) and (b) are plots of the evolution of four-particle eigenstates 
[ipi) and 1^2) respectively, in a field configuration generated according to the 
probability distribution of with X = 0.65 and 9 = 0.267r. 
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Figure 3: Plot of Bio{t), against lattice time t, for X = 0.95, 6 = O.Itt and 
two different one-particle eigenstates for initial states. 
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Figure 4: Plot of log(rc) against log(e) 
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Figure 5: Convergence time vs 9. 
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Figure 6: Plot of log Tc against particle number. 



32 








5000 10000 15000 20000 25000 30000 35000 40000 45000 50000 

Elementary motions 



Figure 7: Plot of C„, against the number of elementary motions, for X = 0.95, 
^ = O.lvr and two different one-particle eigenstates for initial states. The data 
is from the same run as shown in figure 
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Figure 8: Plot of Bio{t), against lat- Figure 9: Plot of Cn, against the num- 
tice time t, for X = 0.95, 9 = O.lvr ber of elementary motions. The data 
and two different one-particle eigen- is from the same run as shown in figure 
states for initial states. |H1 
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Figure 10: Plot of i?io(t), against lat- 
tice time t, for X = 0.95 , 6 = 0.257r 
and two different one-particle eigen- 
states for initial states. 
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Figure 11: Plot of (7„, against the 
number of elementary motions. The 
data is from the same run as shown in 
figure 
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